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We present a 2+1-dimensional lattice model for the copper oxide superconductors and their parent compounds, 
in which both the charge and spin degrees of freedom are treated dynamically. The spin-charge coupling parameter 
is associated to the doping fraction in the cuprates. The model is able to account for the various phases of the 
cuprates and their properties, not only at low and intermediate doping but also for (highly) over-doped compounds. 
We acquire a qualitative understanding of high-T c superconductivity as a Bose-Einstein condensation of bound 
charge pairs. 



1. THE QUESTIONS 

The discovery of the new perovskite supercon- 
ductors [jy opened the path to critical tempera- 
tures (T c ) for superconductivity as high as 140 K, 
surpassing by a factor of seven the highest known 
T c for metals or alloys. In fact the physical mech- 
anism for superconductivity seems to be rather 
different. Traditional superconductors present an 
isotope effect, which calls for a phonon mediated 
mechanism. They also have a gap (of order of 
meV) in the excitation spectrum, as shown by 
the exponential decrease of the specific heat at 
low temperature. All this can be accounted for by 
the BCS theory, where electrons of opposite spin 
and momenta are attractively coupled through a 
phonon exchange. The basic BCS statement is 
the existence of a critical temperature at which 
this attraction simultaneously produces the for- 
mation of bound states of pairs of electrons (with 
binding energy of the order of the gap) and the 
condensation of those pairs into a quantum liquid. 

The new superconductors show an intriguing 
phenomenology : 

• Superconductivity appears in doped, ce- 
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ramie materials, like La2Cu04 doped to 
La2-a;Ba x Cu04, starting at x ~ 0.05. 

• At x=0 it is an insulating antiferromagnet. 

• At large x it may be a normal metal. 

• The doping suppresses the antiferromag- 
netic (AFM) correlation between neighbor- 
ing Cu02 planes. However, large (£ ~ 10a) 
in-plain correlations remain just above T c . 
Also, transport phenomena occur mainly in 
the CuC>2 planes, so everything looks like 
a d=2 problem, with localized spins on the 
Cu and mobile holes on the O ions. 

• Again the superconductivity is charge-2 in 
nature (there is a pairing state). 

• Anomalous normal {i.e. non superconduct- 
ing) state: there are experimental indica- 
tions of a pseudo-gap phase (pair formation 
before quantum liquid condensation). 

All this poses some essential questions [||: 
What is the physical origin of the anomalous nor- 
mal state? How can it be characterized? What 
is the mechanism for high T c superconductivity? 
And what is the pairing state? 
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2. OUR MODEL 

The natural (lattice) field-theoretical approach 
is effective field theory. For instance, the long 
wavelength, low temperature (antiferro)magnetic 
excitations of the the parent compound (undoped 
La2Cu04) are well described by the 0(3) Non- 
Linear cr-model (NLctM) 0. The basic ingredi- 
ents of our proposal are: 

• Treat localized spins (at Cu), and charge- 
carriers (at O) independently and, both, dy- 
namically i.e. not hole moving on a fixed 
spin background. 

• Relativistic lattice field theory, where M.C. 
and Mean Field calculations are feasible. 

• Try to catch the essential features of elec- 
trons strongly coupled to AFM spin back- 
grounds (i.e. not a microscopic model). 

We define a 2+1-dimensional 0(3) lattice 
model, with two flavors of naive fermions: 

s = \ E *M*£w. -*£-/.) (!) 

where <f) = {<j>x,<f>2,<H)\ </> 2 = 1 and ^, * are two- 
component Dirac spinors in d = 2 + 1. The fla- 
vors mimic the two values of the electron spin. 
The matrices <7 M ,Ti are Pauli matrices. The 
spin-fermion coupling y will be related to the dop- 
ing fraction x. In the large-y region, the change 
of variable — * f*,# -> shows that the 
fermionic kinetic term is suppressed by a factor 
1/y. Therefore, at both y — and y — oo we 
recover the NLerM. The insulating undoped ma- 
terial should be represented by the y — > oo limit, 
as suggested by the fermion hopping suppression. 
Thus small x corresponds to large y (a possible 
correspondence is x ~ C 2 /(C 2 + y 2 ), C being a 
constant). 

There is a technical problem for the numeri- 
cal simulation, as <f> are constrained variables not 
belonging to a He group. We solve it by introduc- 
ing three conjugate momenta per spin. We then 
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Figure 1. Phase diagram of the action (g), with 
two fermion families as required by HMO Dashed 
lines indicates the MF calculation, while solid are 
from a MC calculation on a 8 3 lattice. 



obtain the Hamiltonian 

n 

for which equations of motion preserving the con- 
straint and the energy, can be written: 

s s 

4> n = P„ x $„ , P n = -<j> n x ~y-. 

0([>n 

The standard HMC algorithm is now straight- 
forward. The simulation took 16 days of the 32 
Pentium Pro processor parallel computer RTNN 
based in Zaragoza. In fig. |l], we show the phase 
diagram of the model at zero temperature. No- 
tice that it is very similar to the familiar dia- 
gram of (chiral) Yukawa models We name 
their ferrimagnetic phase exotic magnetic (EM), 
to avoid confusion in a condensed matter con- 
text. So far, the undoped material could be rep- 
resented by either of (±KQ,y—oo) in fig. [|. How- 
ever, while doping the material, the represen- 
tative point moves in the plane (k,J/). Doping 
suppresses the AFM order, so the curve n—K(y) 
should start at (—K ,y — oo). In fig. [l], the arrows 
sketch the possible trajectory of a material with 
increasing doping. 
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For the fermionic correlation functions, our re- 
sults are still at the mean-field level. The small- 
y region is in a perturbative regime, with light 
fermions of mass rnt = y{<f>) (a Fermi liquid). In 
the large-y region, the fermionic kinetic term in 
the action is (l/y^^d^ ■ t)W at the MF level. 
Thus, inside the PMS phase ({$ ■ r) = 0), the 
fermions are essentially non-propagating. How- 
ever, the operator effityftyf excites a spin singlet 
of mass rrip air . — 4(y 2 — 3/2), which can get small 
in the PMS phase, that is a light charge-2 spin 
singlet bound state, as in (chiral) Yukawa models 
JHJ. Inside the FM(S) phase, the kinetic term is 
not so strongly suppressed ((</>) ^ 0). Thus, car- 
riers are expected to propagate. In the AFM(S) 
phase the original fermions (and doublers) do 
not propagate. But fermions with wavenumbers 
±7r/2 are found to propagate (strongly resembling 
the Schrieffer pockets) . 




Figure 2. Qualitative temperature (T) and dop- 
ing (x) phase diagram, from the model pi). 



Now, we can qualitatively discuss the behav- 
ior with the temperature (see fig. and ref. 
U). The T=0 axis can directly be read off from 
the arrow line in fig. [|. As we enter the PMS 
phase, we should expect the light, charge-2 bound 
states to be Bose-Einstein condensed, yielding su- 
perconductivity. As the temperature rises, the 



BE condensation will disappear, and the system 
enters a region with heavy single fermions and 
uncondensed light charge-2 pairs. This we ex- 
pect to correspond to the pseudo-gap phase. At 
even higher temperatures, the pairs will break up 
(yielding an insulating phase at the MF level). 
We also obtain Fermi liquid behavior at large x 
(small y). 

3. OUR ANSWERS 

• Physical origin of anomalous normal 
state? A dynamical, antiferromagnetically 
interacting spin background, strongly cou- 
pled to fermions [heavy <f)^> fermions). 

• Characterization of anomalous normal 
state? heavy single charges (0\f r ) and light 
bosonic charge-2 bound states. 

• Mechanism of high-T c ? Bose-Einstein 
condensation of previously formed stable 
pairs. 

• What is the pairing state? A bound 
state of heavy fermions, tied by spin-waves 
in a disordered phase: a PMS-pair. 

Moreover, this pairing mechanism, as well as 
the non-coincidence of pair formation and quan- 
tum liquid condensation, are likely to occur in 
other spin-fermion models. 

REFERENCES 

1. J.G. Bednorz and K.A. Miiller, Z. Phys. B64 
(1986) 184. 



D. Pines, preprint sond-mat /9704102 . 
S. Chakravarty, B.I. Halperin and D.R. Nel- 
son, Phys. Rev. B39 (1989) 2344. 
See, e.g. J. Shigemitsu, Nucl. Phys. (Proc. 
Suppl.) B20 (1991) 515. 
M.A. Stephanov, Phys. Lett. B266 (1991) 
447. 

For further discussions and references, see 
J. L. Alonso, Ph. Boucaud, V. Martin-Mayor, 
A. J. van der Sijs, preprint cond-mat/9706022 
(and preprint SCSC-TR-97-16). 



